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Diel vertical migration (DVM) of zooplankton is a global phenomenon, characteristic of both marine and limnic 
environments. At high latitudes, patterns of DVM have been documented, but rather litde knowledge exists regarding 
which species perform this ecologically important behaviour Also, in the Arctic, the vertically migrating components 
of the zooplankton community are usually regarded as a single sound scattering layer (SSL) performing synchronized 
patterns of migration direcdy controlled by ambient light. Here, we present evidence for hitherto unknown complexity 
of Arctic marine systems, where zooplankton form multiple aggregations through the water column seen via acoustics 
as distinct SSLs. We show that while the initiation of DVM during the autumnal equinox is light mediated, the vertical 
positioning of the migrants during day is linked more to the thermal characteristics of water masses than to irradiance. 
During night, phytoplankton biomass is shown to be the most important factor determining the vertical positioning of 
all migrating taxa. Further, we develop a novel way of representing acoustic data in the form of a Sound Image (SI) 
that enables a direct comparison of the relative importance of each potential scatterer based upon the theoretical con- 
tribution of their backscatter Based on our comparison of locations with contrasting hydrography, we conclude that a 
continued warming of the Arctic is likely to result in more complex ecotones across the Arctic marine system. 

KEYWORDS: diel vertical migration; Arctic; ADCP; acoustics; zooplankton; physical forcing; predator— prey inter- 
actions; light regime; water masses; complexity 
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INTRODUCTION 

Die! vertical migration (DVM) of zooplankton is a char- 
acteristic feature of the world's oceans and lakes, and is 
suggested to be the largest synchronized movement of 
biomass on the planet (Hays, 2003). Since the phenom- 
enon was first described almost two centuries ago, there 
have been numerous studies into both the adaptive sig- 
nificance of this behaviour and its ecosystem conse- 
quences (for two reviews, see Hays, 2003; Ringelberg, 
2010). Research on DVM has focused both on the prox- 
imate and the ultimate explanations of the behaviour (e.g. 
Hays et al., 2001; Lampert, 1989; Ohman, 1990), as well 
as on the implications for the ecosystem regarding 
trophic interactions and the biological carbon pump (e.g. 
Buesselergffl/., 2007). 

In the Arctic, the phenomenon of DVM is rather poorly 
understood, yet it is likely to have an important role to play 
in terms of the fate of carbon sequestration, food web 
interactions and the coupling processes between the 
pelagic and benthic realms, hi general, zooplankton per- 
forms vertical migrations to balance the need to feed close 
to the surface where food is generally most abundant, but 
also to reduce the accompanying risk of being eaten in 
these shallow waters (e.g. Hays, 2003) where illumination 
creates favourable conditions for visual predators (Gliwicz, 
1986). In the high Arctic, with its unique light regime, this 
behaviour has been frequently observed during autumn 
and spring when the day-night cycle is pronounced (e.g. 
Blachowiak-Samolyk et al., 2006; Cottier et al., 2006; 
Falk-Petersen et al., 2008; Wallace et at, 2010), whereas a 
number of studies have failed to find any coordinated verti- 
cal migration during periods of continuous light (e.g. 
Blachowiak-Samolyk et ah, 2006; Fisher and Visbeck, 
1993; Kosobokova, 1978). Recently, Berge et al. (Berge 
et al, 2009, 2012a) demonstrated that active DVM is per- 
formed by multiple species of zooplankton even during 
the polar night at a time of year when the ambient light is 
generally assumed to be insufficient to cue their ascent/ 
descent. Wallace et al (Wallace et al, 2010) identified six 
separate phases of DVM throughout the year, each of 
which had a distinct pattern, suggesting that different 
species were dominating the observed patterns during dif- 
ferent seasons. However, our current state of knowledge 
does not yet permit any reliable identification of the 
species responsible for this behaviour, nor the mechanisms 
regulating the vertical positioning of the organisms. Both 
of these factors are critical for a holistic understanding of 
the patterns and processes governing Arctic marine 
systems, and are major obstacles for any predictive under- 
standing of how ecosystem processes are likely be affected 
by the current warming of the Arctic Ocean and sur- 
rounding shelf seas. 



Most DVM studies from the Arctic (but see, e.g. 
Cisewski et al., 2010; Lavery et al., 2007 for examples in 
other regions of the world) have been carried out using 
either echo sounders (Falk-Petersen and Hopkins, 1981; 
Falk-Petersen et at, 2008; Rabindranath et al, 2011), 
ADCPs (Berge et al, 2009; Cottier et al., 2006; Fisher and 
Visbeck, 1993; Wallace et at, 2010) or zooplankton nets 
(e.g. Blachowiak-Samolyk et al., 2006; Longhurst et al., 
1984), with the robustness of the observations and con- 
clusions limited by the inherent bias and limitations of 
each of the sampling methods. Many of these studies 
have been opportunistic with interpretation of the DVM 
signal being based on incomplete or non-coincident data 
series. From the Arctic, partiy for logistical reasons asso- 
ciated with sampling during the long polar night, most 
studies have adopted an approach involving moored 
ADCPs. There is a long track record of using ADCPs to 
study DVM patterns both in the Arctic and elsewhere 
(e.g. Cisewski et al., 2010; Cottier et al., 2006; Tarling et al., 
2002; Wallace et al, 2010 for some recent examples) with 
the acoustic signal converted to a volume backscatter. 
However, this use of an ADCP is not ideal as the instru- 
ments are generally not calibrated as would be a standard 
echo sounder. However, and despite the concerns regard- 
ing aspects of the data provided by the instruments, they 
are widely used both in moored operations as well as 
fitted onto moving vessels and AUVs (Berge et al., 2012a; 
Schofield et al., 2010). Acoustic instruments moored over 
the entire annual cycle can provide us with a unique tem- 
poral appreciation of DVM patterns. Such a seasonal 
perspective is essential to achieve a more complete and 
quantitative understanding of pelagic community dy- 
namics as well as the role of DVM in zooplankton annual 
routines (Varpe, 2012). 

Although important for our understanding of the 
pelagic system, methods using ADCPs typically emphasize 
the broader migratory picture rather than the detailed 
characteristics of the signal (Wallace et al., 2010). 
Accordingly, and in line with a traditional view of short, 
low-diversity food chains in Polar Regions (Smetacek and 
Nicol, 2005), patterns of DVM in the Arctic are generally 
conceived as a single uniform layer of zooplankton, visible 
in the acoustic data as a Sound Scattering Layer (SSL), mi- 
grating in one synchronized response. While we know 
from other systems that the vertically migrating compo- 
nents of the pelagic community can often be seen as dis- 
tinct SSLs (e.g. Klevjer et al., 2012 from the Red Sea), such 
information and understanding is to a large part missing 
from Arctic systems. Yet, given the importance of current 
changes in primary productivity (Arrigo et al., 2012) and 
sea ice extent (CavaJieri and Parkinson, 2012), there is an 
imminent requirement to understand the processes regu- 
lating vertical migration of zooplankton in the Arctic. This 
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is particularly important if we want to properly account for 
its contribution to carbon cycling (Wallace et ah, 2013) 
through coupling processes between the sympagic, pelagic 
and benthic systems. 

Here, we address the nature of Arctic DVM patterns in 
terms of their complexity, species composition and 
drivers (biotic and abiotic). We present the results from a 
study of DVM in the Arctic collected from four contrast- 
ing sites around the high Arctic archipelago of Svalbard 
(78-80°N) during the autumn equinox, combining net- 
based zooplankton collection with moored acoustic 
instruments at a time of strong diel change in illumin- 
ation. Based on this, we address our two main research 
questions. First, to what extent is there uniformity in the 
patterns of DVM occurring in arctic locations with con- 
trasting hydrography? Second, can variations in DVM 
patterns be related to the physical environment? When 
observing DVM with acoustics, it is inevitable that we are 
also interested in the species composition of any SSL. 
Therefore, net samples were used both to detect organ- 
isms present at different depth layers in the water column 



and as a basis for calculating the theoretical backscatter 
of organisms that de facto are present in the water column. 
Understanding the identity of the migrants provides im- 
portant insights into cues for vertical migration and the 
potential implications of this behaviour for carbon reten- 
tion and export. It is also a vital step towards a more com- 
prehensive understanding of the patterns and processes 
governing Arctic marine systems in an era of strong cli- 
matic change. 



METHOD 

Study area 

The study was carried out from R/V Hehmr Hanssen 
between 16th and 28th of September 2010. The period for 
observations was centred on die equinox to ensure the 
maximum diel cycle in solar irradiance. Data were collected 
at four sites in the Svalbard archipelago: Kongsfjorden, 
Rijpfjorden, BiUeljorden and an Ice Station north of 
Spitsbergen (Fig. 1 and Table I). Kongsfjorden on the 




16°E 



20°E 



24°E 



Fig. 1. Regional map of the study area with stars indicating tlie four sampling sites and labelled Kongsfjorden (KF), Rijpfjorden (RF), Billefiorden 
(BF) and Ice Station (ICE). Bathymetric contours are at 200, 500 and 1000 m with Uie relief enhanced by shading. 
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Table I: (A) Mooring position, depth and deployment time at each site, and the depth layers for MPS 
(multinet plankton sampler) and MIK (Method Isaac Kidd)/WP3 net samples at the four sites (B) 
Acoustic Doppler current profiler (ADCP) mooring arrangement giving the instruments depths, orientation 
and profiling range. 

Mooring deployment 



Site 


Position 


Depth (m) 


and recovery dates/times 


MPS depths (m) 






MIK/WP3 depths (m) 


(A) 
















Kongsfjorden 


78-59.58N 


315 


16.09.1020:00 


300-200, 200-100, 100-50, 50- 


■20, 


20-0 


100-0 m (MIK) 




1 V32.70E 




18.09.1007:00 










Rijpfjorden 


80=1 7. 07N 


212 


19.09.1005:30 


215-100, 100-50, 50-20, 20-0 






50-0 m (MIK) 




22' 16.08E 




20.09.10 16:30 










Billefjorden 


78'39.79N 


186 


27.09.10 06:30 


150-100, 100-50, 50-20, 20-0 






50-0 m (MIK) 




16=45. 42 E 




28.09.10 20:30 










Ice Station 


81 = 10N 


2000 


23.09.10 14:00 


350-200, 200-100, 100-50, 50- 


■20, 


20-0 


50-0 m (WP3) 




11'05E 




25.09.10 02:30 










(B) 
















Site 


Instrumentation 


Depth (m) 


Orientation 


# Data Bins 






Depth range (m) 


Kongsfjorden 


ADCP 1 


96 


Upward 


21 






90-10 




ADCP 2 


191 


Upward 


24 






185-93 




ADCP 3 


194 


Downward 


30 






200-316 




CTD 


31.5 












Rijpfjorden 


ADCP 1 


83.5 


Upward 


18 






77-9 




ADCP 2 


175.5 


Upward 


24 






169.5-77.5 




CTD 


17 












Billefjorden 


ADCP 1 


79 


Upward 


17 






73-9 




ADCP 2 


82 


Downward 


30 






204-88 




CTD 


28.5 












Ice Station 


ADCPl 


6.5 


Downward 


20 






88.5-12.5 




ADCP 2 


162.5 


Upward 


18 






156.5-88.5 




ADCP 3 


164 


Downward 


30 






286-170 




CTD 


31.5 













The CTD was fitted with a PAR sensor. The same instruments and sensors were used at all stations. 



western coast of Spitsbergen is periodically dominated by 
intrusions of warm and saline Atlantic Water (AW) from the 
West Spitsbergen Current (Cottier et ai, 2005; Svendsen 
et ai, 2002) and has been largely free of sea ice during the 
past six winters. It displays many sub-arctic characteristics 
in terms of its water properties and this is reflected in the 
zooplankton species composition (Kwasniewski et ai, 2003; 
Willis et ai, 2006). Rijpfjorden at the northern part of 
Nordaustlandet is predominantly influenced by Arctic 
Water (ArW) (Ambrose et ai, 2006; Soreide et ai, 2010) and, 
in contrast to Kongsfjorden, is covered by sea ice for 6—8 
months of the year. Billefjorden is a silled fjord in the irmer 
part of the Isfjorden complex on the western coast of 
Spitsbergen. It is ice-covered for most of the winter and 
spring each year (NUsen et ai, 2008), and experiences rather 
little influx of AW (Arnkvajrn et al, 2005; NUsen et ai, 
2008). The Ice Station north of Spitsbergen was located 
over deep water (> 2000 m) and within the Marginal Ice 
Zone with a very dense cover of drifting sea ice. In our later 
analysis, water masses are defined according to Cottier et al. 
(Cottier et al., 2005); AW witii T > 3.0°C and S > 34.65, 
ArW with T<1.0°C and 34.30 < S < 34.80, surface 
water (SW) with T > 1.0°C and S < 34.00 and Wmter 
Cooled Water (WCW) having T < - 0.5°C and S > 34.4. 



In order for sampling to coincide with the anticipated 
close coupling between the timing of vertical migrations 
and the timing of the daily cycle of irradiance (Cottier et al., 
2006; TarHng et ai, 2002), all net sampling was conducted 
at a time corresponding to the local sun noon (LSN) (when 
the sun is at its zenith) and the local midnight (LM). LSN 
and LM were determined for each site by their longitude 
with the time of LSN/LM changing by 4 min for each 
degree of longitude. The longitude range in this study 
spans 1 1 degrees (see mooring positions in Table lA) equat- 
ing to a difference in 44 min in the timing of LSN/LM 
between the most westerly (Ice Station) and most easterly 
(Rijpfjorden) sites, despite all sites being within the same 
time zone. To aid comparisons between sites, all data have 
been plotted onto a time axis based on the local time for 
each site, and all the times reported throughout the paper 
are referenced to these local times. 

Acoustics and physical parameters 

At each site, an instrumented mooring was deployed for a 
minimum of 36 h, thereby ensuring that the nets and 
acoustics covered at least one fuU DVM cycle. The moor- 
ings in the shallower fjord sites were anchored to 
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the seabed and retrieved using an acoustic release. At die 
deeper Ice Station, the mooring was suspended from the 
ship, which then followed the path of the ice floe to which 
it was anchored; drifting ~4.5 nautical nules during die 
2-day sampling period. At each site, a Seabird Electronics 
91 1 CTD unit was used to obtain vertical profiles of tem- 
perature and salinity. The unit was fitted with a calibrated 
Seapoint fluorometer to provide Chi a (|xg L profiles as 
a measure of phytoplankton biomass, measured as emitted 
fluorescence at 685 nm. Additionally, a Seabird Electronics 
19+ CTD unit fitted with a Biospherical Instruments 
scalar irradiance sensor, E„ (400-700 nm, photosynthetic 
active radiance, PAR, |jLmol photons m ^ s was placed 
on the moorings. 

Each mooring was fitted with 307.2 kHz RDI Acoustic 
Doppler Current Profilers (ADCP) (Table lA and B). As 
each ADCP is able to sample over a range of about 
100 m, two ADCPs were deployed in Riipfjorden (212 m) 
and BUlefjorden (186 m), and three in Kongsfjorden 
(315 m) and at the Ice Station (2000 m). FuU water column 
coverage was limited by the inability of the ADCPs to 
record data within the blanking distance of 4 m from the 
transducer head and within 1 0 m of the surface. The 
ADCPs were configured to measure the mean echo 
strength from contiguous bins of 4 m and from ensembles 
of 60 acoustic pings at a rate of one ping per second. The 
echo-strength data was then converted to a measure of the 
absolute backscatter (Sv, dB), following the application of 
the SONAR equation outlined in Deines (Deines, 1999, 
but see also Cottier et ai, 2006, and references therein). 
The absolute backscatter over a 36-li period at each site is 
presented in Fig. 2 and Supplementary data. Fig. SI. 
Swimming speeds of the constituent zooplankton are in- 
ferred from the inclination of the SSLs (Cottier et al, 2006). 

We note that the transducers on the ADCPs are not cali- 
brated using standard target techniques (Foote et al, 1987) 
and hence do not provide an accuracy comparable with 
echosounder techniques. Tliis is because backscatter 
strength is not always directly related to the size of a target 
or the abundance of targets, as physiological differences in 
shape, tissue type and orientation will affect the target 
strength of scatterers (e.g. Stanton et al., 1994, 1998a,b and 
references therein). However, careful application of tiie 
methodology following Deines (Deines, 1999) has been 
shown to give good correspondence to echo sounder tech- 
niques (Brierley et al, 2006). Further, organisms that contain 
a swimming bladder such as fish wiU give a much higher 
target strength and scattering value than the zooplankters. 
Accordingly, one would expect that any pelagic fish that 
enters the acoustic beams would give rise to a strong back- 
scatter signal. However, such strong backscatter signals 
from fish are routinely rejected by the ADCP internal data 
quaHty filters (RDI, 1996). 



Zooplankton 

Zooplankton sampling took place close to the fixed moor- 
ings. Mesozooplankton were sampled with Multi Plankton 
Sampler (MPS, Hydro-Bios, Kiel, equipped with five 
closing nets, mesh size = 200 |xm, opening = 0.25 m^). 
The sampling depths were standardized between sites and 
are presented in Table I. To capture any change in the verti- 
cal positioning of zooplankton at the extremes of the did 
cycle of illumination, the MPS sampling was carried out 
during LSN and LM. Macrozooplaiikton (i.e. larger chae- 
tognaths, kriU, ctenophores, ampliipods and pteropods) 
were immediately removed from the sample, measured 
(length, width) and counted, after which the samples were 
stored in borax buffered 4% formaldehyde-in-seawater so- 
lution and transported to land for later analysis. In the lab, 
the smaller remaining size fraction of zooplankton was iden- 
tified to the species level. Sub-samples of 5 mL, obtained 
using a wide-mouthed pipette, were counted until 100 indi- 
viduals of Calams spp. from each sample were identified to 
stage and then measured according to Kwasniewski et al. 
(Kwasniewski et al., 2003). Only stages CIV and above were 
included in analyses, because these stages are by far the 
most dominant at this time of year at the selected sites (e.g. 
Daase et ai, 2008), but also because the contribution to 
acoustic backscatter decreases with organism size. 

Macrozooplankton was collected during one LSN and 
one LM with a MIK net (mesh size =1.5 mm with a filter- 
ing cod end, opening = 3.15 m^) in the upper 50 m (100 m 
in Kongsfjorden). At the Ice Station, with a high concentra- 
tion of sea ice, the larger MIK net was replaced by the 
smaller WP3 net (mesh size = 1000 |xm, opening = 1 m^). 
For all nets, the counts were converted to individuals per 
unit volume (ind. m and are presented in Table II. To 
minimize any influence of spatial heterogeneity of the 
sampled zooplankton community, sequential net hauls were 
taken over a maximum period of 1 h. Length and width of 
up to 100 random individuals in each taxonomic group at 
each site were measured in order to establish a size distribu- 
tion, length-width ratios and standard deviations needed 
for target-strength calculations (see below). 



Identification of Sound Scattering Layers 

For each site, the vertical distribution of acoustic SSLs at 
LSN was identified objectively in the backscatter data 
(Fig. 2). A mean vertical profile of backscatter at LSN was 
calculated from those data spanning a 1 h period; 30 min 
either side of the time of LSN. The SSLs were then 
defined to be the peaks in the backscatter profile where 
values exceeded the overall mean backscatter calculated 
from all data collected during the 36-h measurement 
period. The vertical extent of each SSL was also 
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Table II: Abundance (m of the main 


zooplankton groups 


at LSN and LM and the theoretical 


backscatter (anj 


in brackets calculated for all organisms collected using Eqn (2) and calculated TS values 




Kongsfjorden 




Billefjorden 




axon 


LSN 0-100 m 


LM 0-100 m 


LSN 0-50 m 


LM 0-50 m 


^] 

Mertensia ovum 


0.127 


0.283 


0.994 


0.636 


Krill (small) 


n -1 /I O / 1 AO\ 

U. 1 4o ( — 1 Do) 


1 010/ DD\ 

1 .J 1 2 ( — yy) 


A AAC < 1 00\ 

U.UUb (— 1 zz) 


0 


Krill (large) 


U 


A 

(J 


A A1 Q / 1 no\ 
u.u 1 y I — 1 uz; 


Z.o /o I — oUJ 


Chaetognaths 


1 m -1 / 1 oo\ 

1 .y/ 1 \ — 1 ZZ) 


1 10/1/ 1 no\ 
1 . 1 o4 (— 1 Oo) 


A 00/1 < 1 AO\ 

(j.yy4 (— 1 oy) 


-looni / n"7\ 

1 z.yy 1 ( — y/) 


Themisto abyssorum 


u.uD 1 \ 1 uy/ 


1 .UD4 \ yoj 


A 

u 


A 

u 


Themisto Hbellula 


0 


0 


A ADO < D01 

u.uoo (— yz) 


A OCO f Q 1 1 

u.ybo ( — 0 1 ) 


Other zooplankton 


U.Uoz 


n nno 






Calanus spp. 


"7/1 O / DO\ 


c /I A / oc;\ 


14 ( — 91 ) 


01 Of 001 

z 1 .y ( — by) 


3) 

axon 


Rijpfjorden 




Ice Station 






LSN 0-50 m 


LM0-50m 


LSN 0-50 m 


LM0-50m 


Mertensia ovum 


3.885 


1.032 


0.02 


0 


Krill (large) 


0 


0.356 (-89) 


0 


0 


Krill (small) 


0.134(-109) 


0.153 (-108) 


0.02 (-11 7) 


0 


Chaetognaths 




y.ouo \ yy/ 


n 1 R ( — 1 1 Ri 
u. 1 0 t 110/ 


n R4 f— 1 nQi 
u.c54 \ 1 uy/ 


Themisto libellula 


0.229 (-87) 


0.064 (-93) 


0.18(-88) 


0 


Tiiemisto abyssorum 


0 


0.013(-115) 


0 


0 


Limacina lieiicina 


0.057 


0.013 


0 


0 


Onisimus glaciaiis 


0.006 


0.013 


0 


0 


Other zooplankton 


0.051 


0.102 


0.02 


0.04 


Caianusspp. 


15.5 (-90) 


14.7 (-90) 


3 (-97) 


23.14(-88) 


Paraeucliaeta glaciaiis 


0 


0 


0 


0.28 


Ostracods 


0 


0 


0 


0 


Other zooplankton 


0 


0 


0.02 


0.04 


Apherusa glaciaiis 


0 


0 


0 


0.6 


Metridia sp. 


0 


0 


1.0(-102) 


78.0 (-83) 



(A) Kongsfjorden and Billefjorden, (B) Rijpfjorden and Ice Station. 



determined to be the range in the water column where 
the backscatter at LSN exceeded the overall mean back- 
scatter. We also performed similar analysis for identifying 
SSLs at LM. 



Analyses of acoustic backscatter 

Scattering models have been developed for three main 
types of scatterers (Stanton et al., 1998a,b); fluid-filled 
organisms such as kriU, elastic-shelled organisms such as 
gastropods, and gas-bearing organisms such as siphono- 
phores. To interpret the measured values of absolute back- 
scatter, the target strength (TS, in units of dB) of each 
zooplankton species was calculated based on the randomly 
oriented fluid bent cylinder model (Stanton et al., 1994): 



TS= 10 log <^ 0.087^2^2 /3^i 



l-exp(-87r2/Wc-2) 



cos \ tt/Dc'^ ( " ^ ' + ^■'^ 



where R is the reflection coefficient = 0.058 (Greene et al., 
1998), L the mean body length (m), s the (standard devi- 
ation of I)/L, D the mean body width (m), f^i, = L/D,f 
the acoustic frequency in Hz = 307.2 x 10^, and c the 
speed of sound in water, defined here as 1500 m s The 
theoretical backscatter (Sv, dB) of each taxonomic group 
of zooplankton present in the water column is then esti- 
mated by (see also Cottier et al., 2006; Stanton et al., 
1994): 



Sv(C) = 10%(C) + TS 



(2) 



(1) 



where C is the net-catch estimate of each zooplankton 
group (ind. m ' ). The theoretical volume backscatter at 
abundance level of 1 ind. m ^ was calculated, which is 
numerically identical to the TS and is referred to here as 
Sv(l). Of the species treated in the present study, the 
chaetognaths and the ctenophore Mertensia ovum are most 
likely the two taxa least suited to the bent fluid model 
used for calculating target strength (Stanton et al., 
1998a,b). Neither possesses a hard exoskeleton or cara- 
pace like the other species (all crustaceans), and their 
target strengths are hence most likely overestimated. 
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From the theoretical estimate of backscatter for each 
zooplankton group, we are then able to identify those 
organisms that wiU dominate the measured acoustic 
signals. TS was calculated for each of the most common 
species from the region, with the lower limit in size being 
the older (>CIV) copepodite stages of Calanus spp. at 
around 2.2 mm. For the Ice Station, with its general low 
abundance of zooplankton, we also took Metridia longa 
into consideration due to its relatively high abundance at 
this site. We used the TS value calculated for Calanus spp. 
as an approximation for M. longa as the two taxa resemble 
each other in external morphology and size. For the krill 
species, the size spectrum obtained during the cruise was 
used to sub-divide the population into two classes; large 
(mean length 18 mm) and small (mean length 8 mm) 
krill. All groups examined are hereafter referred to as 
"species". To enable a viable comparison between each 
of the species' individual contributions to the mea- 
sured backscatter levels, both the numerical difference 
(A Sv) in Sv(l) values between the species and the 
linear transformation (according to Stanton et ai, 
1994) of A Sv following A Sv = 10^'^^^"'^ are presented 
in Table III. 

We then calculate a total estimated backscatter for 
each net haul, which we term the Sound Image (SI). The 
SI (using Equation 3) is the sum over all theoretical back- 
scatter values for each species at each site and is based on 
the abundance data from the net hauls presented in 
Table II. 

SI = 10 log(10^'"''° + 10^''2/'° + ■ ■■+10'^''"/'°) dB 

(3) 

where Sv„ is the estimated Sv value (Eq. 2) for species n. 

The SI reflects the relative contribution (in dB) from 
each species to the total estimated backscatter, enabling a 
more direct comparison of backscatter contribution than 
from single species dB values alone. 



The relative contribution to the SI of each species is 
presented as sectors within a pie chart, where the entire 
pie represents the entire SI value. All calculations for the 
SI are based on MIK nets (or WPS at the Ice Station), 
except the values for calanoid copepods which are based 
on MPS samples in the equivalent depth intervals. 

As the MIK net used to sample the macrozooplankton 
could not be closed, the SI value can only be determined 
unequivocally for the upper part of the water column at 
LSN and LM for each site. We make the assumption that 
zooplankton species performing DVM will have different 
abundances in the upper layer during LM compared 
with LSN due to DVM behaviour. Calculations for the 
theoretical Sv values also assume a homogenous distribu- 
tion of scatterers, but by definition an SSL is an aggrega- 
tion of organisms. Therefore, the calculated Sv values 
will under-estimate backscatter values compared with the 
measured backscatter within an SSL. Taking this into 
consideration, and the unavoidable problem of calibra- 
tion (see above and Deines, 1999), the calculated theoret- 
ical Sv values are only comparable within each sample 
and location, and may only be used as an indication of 
relative change in backscatter between LSN and LM. As 
such, we restrict our analyses to the relative variations in 
backscatter values rather than attempting a direct com- 
parison of absolute numbers between the measured and 
the theoretical backscatter values. 



RESULTS 

The time series showing the absolute backscatter (Sv) dis- 
tributions detected by the ADCPs at each site are pre- 
sented in Fig. 2. In each case, the mean Sv value for the 
full deployment was calculated and only those backscat- 
ter values that exceeded the mean value are presented 
using a graduated scale where white represents the mean 



Table III: Calculation of the theoretical target strengths (TS) of the most abundant taxa based on 
measurements of organisms collected across all four stations, the relative difference (A) and its linear 
transformation (linear A). 



Taxon/group 


L 


s 


D 


/3d 


TS 


Sv(l) 


ASv 


Linear A Sv 


Mertensia ovum 


0.018 


0.8889 


0.021 


2.5 


-74.0 


-74.0 


34.6 


2884.0 


Themisto libellula 


0.014 


0.214 


0.003 


4.611 


-80.7 


-80.7 


27.9 


616.6 


Krill (large) 


0.018 


0.111 


0.003 


8.7 


-85.0 


-85.0 


23.6 


229.1 


Themisto abyssorum 


0.006 


0.1667 


0.001 


3.187 


-95.7 


-95.7 


12.9 


19.5 


Krill (small) 


0.008 


0.125 


0.001 


8.8 


-100.0 


-100.0 


8.6 


7.2 


Calanus spp. 


0.00226 


0.24 


0.00083 


2.72 


-102.1 


-102.1 


6.5 


4.5 


Chaetognaths 


0.018 


0.0277 


0.001 


18.7 


-108.6 


-108.6 




1.0 



L: mean body length in m, s: (standard deviation in L]/L, D: mean body width in m, Po. L/D, TS: target strength, calculated based on Eqn (1), Sv (1): 
theoretical backscatter of each taxon using a standard abundance of 1 ind. m"^ expressed as dB, A Sv: numerical difference in Sv (1) between taxa using 
the weakest scatterer as a reference expressed as dB, Linear A Sv: linear transformation of A Sv(1 ) using the equation A Sv = 10'''^''°'. Accordingly, one 
7". libellula has the same scattering potential as 616 Chaetognaths or 137 Calanus spp. 
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Fig. 2. Distribution of the absolute backscatter for values greater than the mean Sv value (in dB) observed over a 36 h period at each of the four 
sampling sites: (A) Kongsfjorden, (B) Rijpfjorden, (C) Billefjorden and (D) Ice Station. The green colour bar ranges from white {the mean Sv value 
over the deployment period) to dark green (the maximum Sv value measure over the deployment period). The full backscatter distribution is given 
in Supplementary data, Fig. SI. LM and LSN are indicated by vertical dashed lines. Ambient solar irradiance {E^j] is overlain in purple and 
normalized to the maximum measured value at each site. The detection of the SSLs at LSN is illustrated in panels to the right of each subfigure. 
The vertical grey line indicates the mean Sv value over the 36 h period, with the Sv profile at LSN overlain (black line). The SSLs are defined as 
those regions of enhanced Sv (relative to the deployment mean value) and depicted by thick grey vertical lines with the grey dot indicating the depth 
of peak Sv in each SSL. 



value. This allows only those regions of high backscatter 
to be shown, which represent the highest abundance of 
scatterers. A normalized measure oi in situ ambient solar 
irradiance, Eq (PAR) from one fixed depth is overlain on 
the backscatter distribution. The right-hand panel in 
each sub-figure shows the mean Sv profile at LSN with 
the SSLs identified as the peaks in backscatter relative to 



the mean Sv value for the deployment. These SSLs are 
numbered for reference, with the circle indicating the 
position of maximum Sv value within each SSL (Fig. 2). 

There is a difference in the number and vertical distribu- 
tion of SSLs between each of the sample sites (Fig. 2). 
Focusing on the distribution of SSLs observed during LSN 
at each site (Table TV), we note the greatest level of 
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Fig. 2 Continued 



complexity in terms of number of SSLs at LSN in 
Kongsfjorden (Fig. 2A). Here, we see four discrete SSLs at 
LSN, the three shallowest being typically 10 m in vertical 
extent while the deepest (SSL4) exceeds 100 m in extent 
but not associated with the seabed. These four SSLs are 
seen to ascend and coalesce during the period of sunset to 
form a single SSL at LM. In Rijpljorden (Fig. 2B), there 
are two very clear SSLs at LSN, one rather diffuse and 
centred around 50 m (SSLl) and a deep one (SSL2) con- 
fined to the bottom 20 m. SSL2 makes a rapid migration 
to the surface during sunset, apparently merging with the 
shallower SSLl to create a single dominant SSL in the 



upper 30 m at midnight. In BUlefiorden, three SSLs can be 
identified at LSN; a rather weak layer at around 50 m 
(SSLl), a well-defined layer extending from around 1 15 to 
150 m (SSL2) and a deep layer confined close to the 
seabed (SSLS). During sunset, all the SSLs are seen to 
ascend giving a strong SSL in the upper 25 m and a more 
dilfuse SSL from the sea bed up to around 1 40 m. Finally, 
the Ice Station was characterized by a generally weak back- 
scatter signal. Nevertheless, we identified two SSLs during 
LSN, one confined to the upper 20 m (SSLl) while the 
second is found below 250 m (SSL2). SSL2 is seen to 
ascend during sunset giving a single SSL at LM. 
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Table IV: Attribution of species to SSLs for each site reporting the depth range, temperature and salinity of 
each SSL. 



Locst ion/SSL 


Dspth r8nQ6 (m) 


Sslin ity 


TsmpsrsturB ( C) 


Asc6nt sp66d (mm s ^) 


SuQcjested miQrators 


Kongsfjorden 












1 


20-40 


34.17 


2.83 


<6 


Mertensia ovum 


2 


60-90 


34.52 


1.27 


<6 


Calanus sp. 


3 


1 20-140 


34.57 


0.94 


<22 


Thsmisto 3iDyssorurnj\^f\\\ 


4 


150-280 


34.78 


2.13 


<22 


Themisto abyssorum/knll 


Rijpfjorden 












1 


30-60 


34.06 


2.72 


<5 


Calanus sp, Chaetognaths 


2 


150-bottom 


34.55 


-0.77 


<22 


Krill 


Billefjorden 












1 


30-50 


33.86 


1.00 


<5 


Calanus sp., Chaetognaths 


2 


120-140 


34.62 


-1.58 


<10 


Krill, Th. libellula 


3 


180-bottom 


34.65 


-1.53 


<20 


Hyperbenthos 


Ice Station 












1 


20-30 


32.77 


-1.38 


<6 


Calanoid copepods 


2 


250- 


34.94 


2.46 


<20 


Themisto libellua 



The swimming speed is based on a simple calculation of the average speed of a SSL as it moves from depth (at LSN) to the surface. 




100, 



Fig. 3. Fluore.scence profiles measured as Chi a (|jLg 1 ) at each site, Kongsfjorden (KF), Rijpfjorden (RF), Billefjorden (BF) and Ice Station (ICE). 
In each profile, the vertical extent of the surface SSL at LM is shown by a dark grey vertical line with a dot indicating the depth of the peak Sv 
within the SSL. 



In Fig. 3, we show the relationship between the con- 
centration of Ghl a (|xg L ') in the upper 100 m and the 
position of the surface SSL at midnight (LM) at each of 
the four sites. The Chi a profiles are comparable at all 
three fjord localities with maximum values measured 
around 20. At each of these sites, the peak of the SSL at 
LM is also located at ^ 20 m, coincident with the peak in 
Chi a. At the Ice Station, the surface SSL also shows peak 
backscatter at 20 m, but the highest Chi a is found imme- 
diately under the ice. In most sites, there appears to be a 
clear association between the vertical position of the 
zooplankton and the location of the Chi a maximum 
(Chi 



The peak backscatter in the surface at LM is located 
precisely at the depth of maximum density gradient, the 
pycnocline (Fig. 4) with the SSL distributed through the 
lower part of the pycnocline at all sites. During daytime 
(LSN), the shallowest scattering layer (SSLl at all sites) 
generally also occupied the lower part of the pycnocline, 
the Ice Station being the exception. Deeper SSLs were 
not associated with any significant density gradients. 
Rather, we see from Fig. 5 that LSN SSLs were located at 
inflections in the temperature -salinity curves. By defin- 
ition, the inflection points correspond to layers in the 
water column bounded top and bottom by gradients in 
temperature and/or salinity. Consequendy, an SSL 
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Fig. 4. Profiles of the density gradient at each site, Kongsfjorden (KF), Rijpfjorden (RF), Billefjorden (BF) and Ice Station (ICE). SSLs identified 
during LSN (light grey vertical lines with dots indicating the depth of peak Sv) and LM (dark grey lines with dots at peak .Sv). 



located at an inflection point is within a layer of water 
that has well-defined temperature and salinity character- 
istics and which are distinct from any layers above and 
below. Water masses represent the broad definitions of 
these characteristic temperature and salinity relation- 
ships. Kongsfjorden shows the greatest complexity in its 
hydrographic structure (Fig. 5A), with multiple inflec- 
tions in the temperature— salinity relationship indicative 
of distinct layers of water. The deep SSL4 is within the 
Transformed A\y being the warmest and most saline 
water in the fjord, while the intermediate depth SSL2 
and SSLS are located within interleaving layers of 
water. In Rijpfjorden, Fig. 5B, the deep SSL2 is in the 
coldest water mass of WCW which is a product of sea 
ice formation the previous winter (Nilsen et al., 2008). It 
is similar in Billefjorden, Fig. 5B, where the intermedi- 
ate and deep SSLs are both within the WCW and there- 
fore not separated by any gradients in the hydrography. 
At the Ice Station, the temperature -salinity distribution 
shows a classic form for the arctic halocline with cold 
and fresh SWs overlying the deeper, saline Atlantic 
layer (Rudels et al., 1996, 2005). We see that the deep 
SSL2 is located precisely at the warmest and most 
saHne water mass (TAW) after a descent of over 250 m 
from the surface. 



The zooplankton community 

The zooplankton community in the upper 50 m (100 m 
in Kongsfjorden) was sampled at all four sites during 
both LSN and LM (catch data summarized in Table II). 
Calanoid copepods were numerically the most dominant 
group at all sites, with Calanus spp. dominating within the 
fjords and Metridia longa at the Ice Station. The second 
most dominant group (numerically) was chaetognaths, 
followed by Mertensia ovum, kriU and Themisto spp. 
(Table II). Themisto libellula were observed by divers 
(J. Berge, pers. Obs.) to be present at the Ice Station 
under the sea ice during daytime. 

Lowest backscatter values were estimated for the chae- 
tognaths and Calanus spp., highest for the ctenophore 
Mertensia ovum. The ratios of echo-energies (A Sv, dB) and 
their linear transformations are presented for all taxa 
examined in Table III. In terms of theoretical backscatter 
and the resulting SI in Fig. 6, Calanus spp. and TT libellula 
were the two most dominant taxa at all the sites, with the 
exception of kriU during LM in Rijpfjorden and 
Billefjorden and Metridia longa during LM at the Ice 
Station. In Kongsfjorden (Fig. 6A), Calanus spp. dominated 
the SI both during LSN and LM, with T abyssorum increas- 
ing in relative importance during LM. In both Rijpfjorden 
and Billefjorden (Fig. 6B and C), a comparable pattern of 
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Fig. 5. Temperature— salinity distributions in (A) Kongsfjorden (thick solid Kne) and Rijpfjorden (dashed line), (B) BiUeljorden (dotted line) and Ice 
Station (thin black line). For each site, the occurrence of the SSLs at LSN is shown by the hght grey portion of each line and with the dot indicating 
the depth of peak Sv within each SSL. Characterization of water masses follows Cottier et al. (Cottier ei ai, 2005) and Nilsen et al. (Nilsen et ai, 
2008): S\y Surface Water; UN, Local Water; WC\y Wmter Cooled Water; W, Intermediate Water; TAW; Transformed Atlantic Water; /W, Atlantic 
water The left-hand figure shows the full temperature— salinity distribution while the right-hand subfigure highhghts tliat part of the profile falling 
within the I\\; TAW' and AVV^ domains. Thin grey lines in the left-hand figure are isopycnals ranging from 102.5 to 1028 kg m ^ 



Calanus spp. and TT libellula dominating the SI during LSN 
occurred. At both these sites, the relative contribution of 
krill increased during LM, indicative of a strong krill 
DVM. At the Ice Station, where very few organisms were 
collected in the net samples, the SI during LSN was domi- 
nated by TT libellula, whereas the SI during LM was totally 
dominated by the calanoids {M. longa and Calanus spp.). 
Importandy, for all calculations of SI, the contribution by 
Mertensia ovum was disregarded, as this species is poorly 
suited to the bent fluid model used for calculating target 
strength (Stanton et ai, 1998a,b) and would therefore dom- 
inate the calculation of theoretical backscatter. 

Diel change in abundance of each taxon suggests that in 
Kongsfjorden (Fig. 6A, right panel) kriU, Mertensia ovum and 
T ahyssomm were performing DVM in the upper 100 m. In 
Rijpfjorden (Fig. 6B, right panel), krill and chaetognaths 
show a diel displacement in abundance that is consistent 



with a DVM behaviour In BiUefjorden (Fig. 6C, right 
panel), all taxa except Mertensia ovum had a higher abun- 
dance in the upper 50 m during LM compared with LSN. 
At the Ice Station (Fig. 6D, right panel), Metridia longa, 
Calanus spp. and chaetognaths all had a higher abundance 
in the upper 50 m during LM than LSN. 



DISCUSSION 

Our current understanding of DVM in the Arctic during 
day- night cycles as observed by acoustic techniques is 
characterized by a rather simple ascent/ descent pattern 
of a single aggregation of zooplankton, or SSL (e.g. Berge 
et al., 2009; Cottier et ai, 2006; WaUace et ai, 2010). The 
distribution of SSLs recorded at LSN (Fig. 2) indicates that 
at all the locations in this study zooplankton migrate down 
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Fig, 6. Sound Images and relative zooplankton abundance at LSN and LM for (A) Kongsfjorden (0— 100 m), (B) Rijpljorden (0-50 m), (C) 
Billeljorden (0—50 m) and (D) Ice Station (0-50 m). The Sound Images are represented as pie charts showing the relative contribution of each 
species to the theoretical backscatter value determined from catch data with the theoretical backscatter value given. Based on the likely 
overestimated TS value oiMertensia ovum, this species is not part of the SI calculation (see Discussion for details). The relative abundance of scatters 
at LSN (white bars) and LM (black bars) in the surface are given as a bar chart (right panels). 



to depth during daytime with a subsequent migration into 
the Chi drmoL layer during the night. This pattern is consist- 
ent with classical ideas on DVM (see Hays, 2003 for a 
review). Also, all vertical migration events appear to be 
synchronized in time, which is a strong indication of a 
light mediated response. Nevertheless, the acoustic back- 
scatter data from this study challenge the classical concept 
of simple DVM comprising a single migrating SSL. 
Instead, our data show a more complex DVM process 
with the formation and migration of multiple layers of zoo- 
plankton distributed through the water column at each 
study site. This evidence for complex patterns of DVM is 
the core result of this study. 



Clil a, with the peak of backscatter typically coincident 
with the Chi flmax- Colocation between the shallow SSL 
and the increased availability of food is a strong indica- 
tion of the ascent of zooplankton in the DVM cycle being 
a biologically driven process (e.g. Ringelberg, 2010). The 
SSL and Chi Omax are found within the pycnocline 
(Fig. 4), a phenomenon that has been previously observed 
(Basedow et ai, 2010; Dekshenieks et al, 2001; HoUiday 
et al, 2010) and used for modelling the vertical migration 
of copepods (Wallace et al, 2013). Gonsequendy, we inter- 
pret the consistent observation of surface SSLs at night as 
being zooplankton targeting the enhanced concentration 
of phytoplankton within the pycnocline. 



Sound Scattering Layers at LM (Fig. 3) 

Our observations show that at all sites, there is a dominat- 
ing SSL near the surface at LM, supporting the classical 
view of DVM behaviour Importantly, the shallow SSL is 
formed from a coalescing of deeper SSLs that ascend to 
the surface at dusk. In Fig. 3, we show that the depth of 
the surface SSL is coincident with the region of enhanced 



Sound Scattering Layers at LSN 
(Table IV and Fig. 4) 

Shallow 

The shallowest SSL at LSN was found just below the 
pycnocline having descended from the region of Chi flmax 
(Fig. 4). It is likely that chaetognaths were a key constituent 
of the surface SSL at those sites where they were found in 
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high abundances, e.g. in Rijpfjorden and BUlefjorden. 
DVM by chaetognaths is a known behaviour (Pearre, 
1973; Sweatt, 1980) and the swimming speeds (Table IV) 
are appropriate. A candidate species in the surface SSL at 
Kongsfjorden was Ad. ovum which was generally found 
shallower than 1 00 m at similar densities during both LSN 
and LM. Moreover, it is reported to be restricted to the 
upper surface layer during autumn (e.g. Raskolf et al., 
2005). At the Ice Station, the surface SSL was found pre- 
cisely within the pycnocKne and most likely dominated by 
T libellula (SI, Fig. 6D). Interestingly, however, and based in 
particular on the net samples (Table IIB), the shallow layer 
during LM is composed also of Metridia longa and the 
ice-associated Apherusa glacialis (Berge et al., 2012b). Thus, 
the constant and slightly higher backscatter signal observed 
at the surface layer may have been the presence of ice 
fauna, moving out of and into the ice at diflFerent times of 
the day. Apherusa glacialis was found in nets only during night 
time, indicating movement away from the ice and into the 
pelagic zone under cover of darkness (Berge et al., 2012b). 
Deep 

Our data show that density gradients at the depth of the 
deepest SSLs are negligible. Rather, the depth of the deep 
SSLs at LSN shows a link to the water properties, particu- 
larly temperature. We note that in both BUlefjorden and 
Rijpfiorden, the deepest SSLs migrated into the coldest 
water, WC\y while at the Ice Station and Kongsfjorden 
the deep SSL is located within the thermal maximum of 
AWs (Fig. 5A and B). We conclude that the vertical pos- 
ition of the SSLs represents a balance between predation 
risk, thermal optima and distance to the food-rich layers. 

In terms of species composition (Table IV) of the deep 
SSLs, the relatively high ascent rates of the deep migrat- 
ing SSLs suggest that kriU and amphipods may be the 
main contributors. In both Rijpfjorden and Billefjorden, 
both groups exhibit a clear displacement in their abun- 
dance during LSN (Fig. 6B and C). Indeed, Themisto spp. 
and kriU have previously been shown to perform DVM 
(e.g. Fortier et al., 2001; Kaartvedt, 2010) over long verti- 
cal distances (Berge and Nahrgang, 2013), and krill and 
Themisto spp. distributions and their contributions to the 
SI match the deep migrating SSLs at all four localities. 
We also saw evidence for deep SSLs undergoing partial 
ascent at midnight in both Rijpfjorden and Billefjorden 
(Fig 2B and G). 

The Sis for LSN and LM in Kongsfjorden are rather 
similar and differ only in the additional contribution 
from Themisto ahyssomm at LM (Fig. 6A). Accordingly, one 
of the deep migrating layers in Kongsfjorden is most 
likely formed by the amphipod T abyssorum. Of the taxa 
collected by the nets at the Ice Station, the chaetognaths 
and calanoid copepods both have vertical distributions in 
accordance with the migrating deep SSL (Table IIB). 



However, the abundances of calanoid copepods were 
high and they dominated the SI during LM (Fig. 6D). 
Further, the reported swimming speeds of chaetognaths 
(Hardy and Bainbridge, 1954) do not correspond with the 
ascent velocities observed for the SSL while the high swim- 
ming speed previously reported for Metridia spp. (up to 
25 mm s Miller, 2004) is within the observed range SSL 
(Table IV). Importantly, however, the lower limit of the 
ADCP in our study was around 300 m, and the SSL prob- 
ably started its ascent from below this depth. Hence, con- 
clusions concerning the ascent rate and absolute depth 
distributions of the SSL during LSN are incomplete. 
Intermediate 

Kongsfjorden and Billefjorden both show SSLs during 
LSN at intermediate depths (Figs 2 and 4). In 
Kongsfjorden, the surface SSLl and the intermediate 
SSL2 could both be assigned to calanoid copepods 
which showed high abundance and theoretical backscat- 
ter within the upper 1 00 m during LSN and LM. 
Further, given that Calanus spp. dominated the SI at both 
LM and LSN (Fig. 6A), indicative of migration being 
restricted to the upper layers, and the reported swimming 
speed of C. finmarchicus of 4.2-18.3 mm s ' (Mauchline, 
1998), the intermediate SSL2 in Kongsfjorden (Fig. 4, 
Table IV) is likely to be composed of calanoid copepods. 
Calanoid copepods have been shown to constitute up to 
80% of the total mesozooplankton biomass in Arctic 
region (Soreide et al., 2008), and they have also been sug- 
gested as candidates for acoustically detected DVM pat- 
terns from two of the four study sites used herein (Berge 
et al., 2009; Cottier et al., 2006). In Billefjorden (Grigor 
et al., 2014), it is likely that the intermediate SSL2 is com- 
posed of taxa migrating from the surface (krill and T libel- 
lula. Table IV) being distinct from the deep SSL3 
associated with deep migrating hyperbenthic taxa. 
Despite the rigorous sampling effort around both LSN 
and LM, the migratory species within these SSLs cannot 
be identified unequivocally as the sampling protocols 
were designed with the expectation of single rather than 
multiple SSLs. 

Environmental factors driving complex 
DVM 

DVM is a ubiquitous feature of mid-latitude oceans 
whereby the zooplankton population resides at depth 
during daylight hours and moves en masse towards the 
surface during darkness (e.g. Hays, 2003; Fortier et al, 
2001). While the proximate cue for such behaviour is 
widely accepted to be changes in illumination, specifical- 
ly the absolute and relative light intensity (e.g. Fortier 
et al, 2001; Ringelberg and Van Gool, 2003; TarUng 
et al., 2002), the vertical extent of migration is often 
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interpreted in terms of biological processes; mainly relat- 
ing to community composition, predation and feeding 
(Fortier et ai, 2001; Ringelberg, 2010; Tarling et al., 
2014). Our result concur with these general features in 
that we observe a strong vertical displacement of biomass 
that is correlated with the diel change in solar illumin- 
ation across all four examined stations. Importandy 
however, although the ascent/ descent may be initiated 
by light, the final "resting" depth appears to be deter- 
mined by hydrography rather than light. We suggest that 
this is the result of a trade-off between predator avoidance 
and a thermal optimum for particular species. Therefore, 
a site with rather complex hydrography is likely to see a 
more complex DVM pattern. This is supported by the 
fact that the number of SSLs at each site varies, while the 
zooplankton community composition remains quite 
similar. Indeed, the site with the fewest number of differ- 
ent scatterers (Table II, Fig. 6), Kongsfjorden, is also the 
site with the most complex hydrography and shows the 
highest number of SSLs. 

Other investigators have noted the relationship 
between zooplankton community structure and the local 
hydrography (e.g. Schulz et al., 2012; Willis et al., 2006, 
2011 for some recent examples). Lawson et al. (Lawson 
et al., 2004) related the vertical posiuoning during DVM 
of macrozooplankton (krill) to the occurrence of distinct 
water masses. However, that study regarded the vertical 
migrating community as a single entity, not a complex 
community segregated into several SSLs that might 
display affinities towards different water masses (see also 
Ashjian et al., 2008). Laboratory experiments investigating 
DVM in estuarine systems have shown that zooplankton 
will adjust their vertical position or migratory behaviour 
in response to the presence of haloclines (Lougee et al., 
2002). The explanation for this has been proposed to be 
one of energetics, with less energy required to remain 
suspended in more saline dense water leading to accumu- 
lations of organisms in and below the halocline (Lougee 
et al., 2002). Similar separation of species has been 
observed in relation to the thermocline, with C. finmarchi- 
cus occupying colder water below and C. helgolandicus in 
the warmer fresher water above a thermocline (Williams, 
1985). Other experimental investigations have deter- 
mined that the vertical distribution of zooplankton is in 
response to density rather than salinity (Harder, 1968). 

Our observations (Figs 4 and 5) challenge the import- 
ance of density as a driver of vertical structure (Harder, 
1968), with multiple SSLs located in regions of uniform 
density in at least two of the fjords (Kongsfjorden and 
Billefjorden, Fig. 4). However, our data show that the 
hydrographic structure of the water column does appear 
to correlate with both the number and the position of 
SSLs. In Kongsfjorden, with the highest number of SSLs, 



the three deepest SSLs are associated with distinct tem- 
perature and salinity characteristics (Fig. 5A). Precise ver- 
tical positioning is seen in Rijpfjorden where the deep 
SSL matches the upper boundary of the WCW layer, and 
at the Ice Station where the deep SSL matches the 
thermal maximum of the AW layer (Fig. 5A and B, re- 
spectively). Another indication for a physical forcing of 
the DVM patterns comes from a comparison between 
the data from Kongsfjorden in 2002 (Cottier et al., 2006) 
and 2010 (data herein). Both data sets were collected 
using ADCPs of the same frequency at the same site and 
time of year While the data from 2002 showed only one 
vertically migrating SSL, the number had increased to 
four in 2010. During this transition period, the macro- 
zooplankton species present in the fjord remained rela- 
tively stable (Hop et al., 2002; Kwasniewski et al., 2012), 
although it should be noted that there has been a con- 
spicuous shift in the relative abundance and dominance 
of these species (Kwasniewski et al., 2012). The hydrog- 
raphy, however, has changed dramatically since 2002, 
with pronounced and episodic influx of AW as well as 
ice-free conditions since 2006 (Berge et al., 2005; Willis 
et al., 2008; Wallace et al., 2010, pers. Obs.), supporting 
the hypothesis that a physical forcing of DVM patterns is 
more important than previously assumed (e.g. Klevjer 
etal., 2012). 

A potentially important environmental factor behind 
the observed DVM patterns is the differences in predator 
distributions between the study sites; as also suggested in a 
comparative study of DVM patterns in Kongsfjorden and 
Rijpfjorden (Wallace et al., 2010). A number of studies 
have shown that the presence of specific predators can 
cause rapid behavioural responses among prey (BoUens 
and Frost, 1989a,b,1991; Fortier et at, 2001; TarKng et al., 
2002). Important predators such as the polar cod (Boregadus 
saida) or little auks {Alle alle) are known to be common at all 
of the study sites (FaUc-Petersen et al, 2008; Renaud et ai, 
2012), and might be an important mediator of the distribu- 
tion of zooplankton. Other and more Atlantic species such 
as the Atlantic cod {Gadus morhud) or capelin {Malktus villo- 
siis) are more dominant in Kongsfjorden compared with 
the other study sites (Berge et ai, in press; Nahrgang et at, 
2014; Renaud et ai, 2012). Although not observed in this 
study, Oilman (Ohman, 1990) showed a reversed DVM of 
herbivores in the presence of both tactile and visual preda- 
tors. However, if visually searching predators such as the 
polar cod or capelin exert a high predation risk, one would 
expect size-dependent migrations in the prey, with larger 
individuals (most easily spotted by the predator) migrating 
to greater depth than smaller prey (Hays, 2003). Such a 
size-dependent migration is partiy seen in our study, with 
the larger macrozooplankton constituting the deeper 
SSLs, although it is impossible to isolate a potential 
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predator avoidance mechanism from merely a swimming- 
speed dependent segregation of migrating zooplankton. 
Other potential predatory species known to be important 
in the system include Themisto spp. (Kraft et ai, 2013, this 
study) and chaetognaths (FaDcenhaug, 1991) which are 
known to feed on calanoid copepods (Falkenhaug, 1991). 

Characteristics of SSLs and conclusion 

We present evidence for the occurrence of distincdy dilfer- 
ent patterns of DVM between sites that are relatively close 
geographically with a zooplankton community comprising 
the same set of species (Fig. 6), but with contrasting hydro- 
graphic characteristics. Calanus spp. has previously been 
identified as a major contributor to the vertically migrating 
component of the zooplankton community during late 
summer and autumn (Blachowiak-Samolyk et ai, 2006; 
Cottier et ai, 2006; Falk-Petersen et al., 2008). We, however, 
based upon our comparison between acoustic and net 
data conclude that the acoustic backscatter from calanoid 
copepods is typically overwhelmed by the signal from the 
larger zooplankton species, kriU and Themisto spp. in par- 
ticular, when both are present in the water column. Also, 
we argue that the differences in the hydrographic struc- 
tures, particularly the local temperature characteristics can 
give rise to complexity in DVM patterns. Based on the dis- 
tribution of SSLs in Kongsfjorden and the Ice Station, 
organisms do not appear to be migrating down to colder 
water in which a lowered metabolism might optimize an 
energy budget (Ikeda, 1985). The deepest migrating SSLs 
at both of these stations migrate down through colder 
water masses to temperature maxima. Similarly, the verti- 
cal positioning of the SSLs seems only weakly correlated 
with salinity, although all of the deeper SSLs are associated 
with salinity levels above 34.50 PSU. 

Our study has revealed a complexity of polar marine 
systems that has been poorly acknowledged. First, we 
have documented the existence of multiple SSL in 
regions previously believed to be characterized by simple, 
uniform patterns of DVM (e.g. Berge et al., 2009; Cottier 
et al., 2006; Wallace et al., 2010). Second, we document 
the distincdy different patterns observed across the sites 
examined (Figs 2 — 5), indicating that autumn DVM in 
the Arctic comprises a complex behaviour by a wide 
variety of planktonic taxa. The different taxa performing 
DVM are likely to have distinct feeding preferences, 
migration-depth patterns and faecal pellet properties. 

As there appears to be a positive correlation between 
complexity of hydrography and number of SSLs, any 
long-term change in Arctic hydrography may effectively 
lead to marked changes in both behaviours of individual 
species as well as in overall ecosystem function. There are 
examples of on-going changes linked with e.g. the 



amount of inflowing /W, alterations of freshwater input 
and local mixing of water masses, potentially leading to 
changes in pelagic complexity and community segrega- 
tion as different scatterers position themselves at different 
depths throughout the water column. Ultimately, this sug- 
gests that while DVM of Arctic zooplankton could have 
significant implications on retention and export of 
organic and inorganic compounds through the biological 
carbon pump on Arctic shelves, a more thorough insight 
into which species are performing the DVM is needed in 
order to fully understand their impact. 

SUPPLEMENTARY DATA 

Supplementary data can be found online at http://plankt 
.oxfordjournals.org. 
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